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Large order behavior in perturbation theory of 
the pole mass and the singlet static potential 
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Abstract. We discuss upon recent progress in our knowledge of the large order behavior in pertur- 
bation theory of the pole mass and the singlet static potential. We also discuss about the renormalon 
subtracted scheme, a matching scheme between QCD and any effective field theory with heavy 
quarks where, besides the usual perturbative matching, the first renormalon in the Borel plane of the 
pole mass is subtracted. 



<N : MASS NORMALIZATION CONSTANT 

On 

r^- ■ In this paper, we review some results obtained in Ref. [1]. 

The on-shell (OS) or pole mass can be related to the MS renormalized mass by the 

series 

oo 

/«os = m m + Y, r n a" + ' , ( 1 ) 

Ph| where the normalization point v = my^ is understood for m-y^ and the first three 

coefficients tq, r\ and r 2 are known [2] (a? = of (v), where rii is the number of light 
fermions). The pole mass is also known to be IR finite and scheme-independent at any 
finite order in a s [3]. We then define the Borel transform 



dt e- t/a ° B [m OS ] (0 , B [m 0S ] (t) = £ r n — . (2) 

o " =0 

The behavior of the perturbative expansion of Eq. (1) at large orders is dictated by the 
closest singularity to the origin of its Borel transform, which happens to be located at 
t = 27T./Po, where we define 



v- 



da s 
dv 



- 2 «>S + m) 2+ "l 



Being more precise, the behavior of the Borel transform near the closest singularity at 
the origin reads (we define u = jjg£) 

B[m os ] (t(u)) = N m \ - +b (1 +ci(l -2u) +c 2 (l -2u) 2 + ■ ■ •) + (analytic term), 

(3) 
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n- 



where by analytic term, we mean a piece that we expect it to be analytic up to the next 
renormalon (u = 1). This dictates the behavior of the perturbative expansion at large 
orders to be 

. /poVr(n+i+fc) A b b(b-\) \ ^ 

The different b, c\, C2, etc ... can be obtained from the procedure used in [4] (see [4, 1] 
for the explicit expressions). We then use the idea of [5] and define the new function 

oo 

D m (u) = £DL"V = (l-2 M ) 1+fo fl[mos](^)) (5) 

= N m V (l + c\ (1 - 2u) + c 2 (l - 2uf + •••) + (1 - 2u) 1+fo (analytic term) . 

This function is singular but bounded at the first IR renormalon. Therefore, we can 
expect to obtain an approximate determination of N m if we know the first coefficients of 
the series in u and by using 

N m v = D m (u= 1/2). (6) 

The first three coefficients: D m , D m and D m are known in our case. In order the 
calculation to make sense, we choose v ~ m. For the specific choice v = m, we obtain 
(up to 0(u 3 ) with u - 1/2) 

N m = 0.424413 + 0.137858 + 0.0127029 = 0.574974 (n f = 3) (7) 

= 0.424413 + 0.127505 + 0.000360952 = 0.552279 (n/ = 4) 
= 0.424413 + 0.119930-0.0207998 = 0.523543 (n/ = 5). 

The convergence is surprisingly good. The scale dependence is also quite mild (see [1]). 

By using Eq. (4), we can now go backwards and give some estimates for the r n . 
They are displayed in Table 1. We can see that they go closer to the exact values of r n 
when increasing n. This makes us feel confident that we are near the asymptotic regime 
dominated by the first IR renormalon and that for higher n our predictions will become 
an accurate estimate of the exact values. In fact, they are quite compatible with the results 
obtained by other methods like the large |3o approximation (see Table 1). 

We can now try to see how the large |3o approximation works in the determination of 
N m . In order to do so, we study the one chain approximation from which we obtain the 
value [6] 

A4 largepo) = ^ei= 0.976564. (8) 

By comparing with Eq. (7), we can see that it does not provide an accurate determination 
of N m . This may seem to be in contradiction with the accurate values that the large |3o 
approximation provides for the r n (starting at n = 2) in Table 1 . Lacking of any physical 
explanation for this fact, it may just be considered to be a numerical accident. In fact, the 
agreement between our determination and the large (3o results does not hold at very high 
orders in the perturbative expansion, whereas we believe, on physical grounds since our 
approach incorporates the exact nature of the renormalon, that our determination should 
go closer to the exact result at high orders in perturbation theory. Nevertheless, the large 
Po approximation remains accurate up to relative high orders. 
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TABLE 1. Values of r n for v = m^. Either the exact result, the estimate 
using Eq. (4), or the estimate using the large po approximation [7]. 



fn — r n/m MS 


h 


h 


h 


h 


h 


exact (iif = 3) 


0.424413 


1.04556 


3.75086 


— 


— 


Eq. (4) (n f = 3) 


0.617148 


0.977493 


3.76832 


18.6697 


118.441 


large p («/ = 3) 


0.424413 


1.42442 


3.83641 


17.1286 


97.5872 


exact (nf = 4) 


0.424413 


0.940051 


3.03854 


— 


— 


Eq.(4)(n/ = 4) 


0.645181 


0.848362 


3.03913 


13.8151 


80.5776 


large Po («/ = 4) 


0.424413 


1.31891 


3.28911 


13.5972 


71.7295 


exact (nf = 5) 


0.424413 


0.834538 


2.36832 


— 


— 


Eq.(4)(n/ = 5) 


0.706913 


0.713994 


2.36440 


9.73117 


51.5952 


large Po («/ = 5) 


0.424413 


1.21339 


2.78390 


10.5880 


51.3865 



STATIC SINGLET POTENTIAL NORMALIZATION CONSTANT 

One can think of playing the same game with the singlet static potential in the situation 
where Aqcd "C 1/r. Its perturbative expansion reads 



n=0 



(9) 



The first three coefficients r , V , , and V_ 2 are known [8]. At higher orders in 
perturbation theory the log dependence on the IR cutoff \ us appears [9]. Nevertheless, 
these logs are not associated to the first IR renormalon (see [1]), so we will not consider 
them further in this section. We now use the observation that the first IR renormalon of 
the singlet static potential cancels with the renormalon of (twice) the pole mass. We can 
then read the asymptotic behavior of the static potential from the one of the pole mass 
and work analogously to the previous section. We define the Borel transform 

oo 

dte-^B[vi 0) ](t) , B[vi%) ee £ Vjf- (10) 

n=0 n - 

The closest singularity to the origen is located at t = 27r/J3o- This dictates the behavior 
of the perturbative expansion at large orders to be 



^**(i)'^K 



-ci + - 



b(b-V 



: + b) (n + b)(n + b-l 
and the Borel transform near the singularity reads 



-C2 + - 



(ID 



B[V s { \t(u))=N v v 



1 



(l-2u) l+b 



(1 + c\ ( 1 - 2u) + c 2 ( 1 - 2m) 2 H ) + (analytic term) . 

(12) 
In this case, by analytic term, we mean an analytic function up to the next IR renormalon 

at u — 3/2. 
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0.345222 = 


-1.10661 


(n/ = 


= 3) 


0.329356 = 


-1.07729 


(n/ = 


= 4) 


0.295238 = 


-1.04175 


(n/ = 


= 5) 



As in the previous section, we define the new function 

D v (u) = £l#V=(l-2») 1+ *fl[V 5 (0) ] (*(>)) (13) 

n=0 

= 7Y V V (l+ci(l - 2«) +c 2 (1 - 2w) 2 + •••) + (1 - 2m) 1+fo (analytic term) 

and try to obtain an approximate determination of Ny by using the first three (known) 
coefficients of this series. By a discussion analogous to the one in the previous section, 
we fix v = 1/r. We obtain (up to 0(u 3 ) with u = 1/2) 

N v = -1.33333+0.571943-0.345222= -1.10661 (n/ = 3) (14) 

= -1.33333+0.585401 
= -1.33333+0.586817 

The convergence is not as good as in the previous section. Nevertheless, it is quite 
acceptable and, in this case, apparently, we have a sign alternating series. In fact, the 
scale dependence is quite mild (see [1]). Overall, up to small differences, the same 
picture than for N m applies. 

So far we have not made use of the fact that 2N m + Ny = 0. We use this equality as 
a check of the reliability of our calculation. We can see that the cancellation is quite 
dramatic. We obtain 

2 2N m +N v _ 
2N m -N v 

We can now obtain estimates for V}^ by using Eq. (11). They are displayed in Table 
2. Note that in Table 2 no input from the static potential has been used since even Ny 
have been fixed by using the equality 2N m = —Ny. We can see that the exact results are 
reproduced fairly well (the same discussion than for the r n determination applies). This 
makes us feel confident that we are near the asymptotic regime dominated by the first IR 
renormalon and that for higher n our predictions will become an accurate estimate of the 
exact results. The comparison with the values obtained with the large Po approximation 
would go (roughly) along the same lines than for the mass case, although the large |3o 
results seem to be less accurate in this case (see Table 2). 

In order to avoid large corrections from terms depending on v us , the predictions should 
be understood with v us = 1/r. 

RENORMALON SUBTRACTED SCHEME 

In effective theories with heavy quarks, the inverse of the heavy quark mass becomes 
one of the expansion parameters (and matching coefficients). A natural choice in the 
past (within the infinitely many possible definitions of the mass) has been the pole mass 
because it is the natural definition in OS processes where the particles finally measured in 
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0.038 


,n/ = 


= 3 


0.025 


,n/ = 


= 4 


0.005 


,n f = 


= 5 



TABLE 2. Values of V s \„ with v = l/r. Either the exact result (when avail- 
able), the estimate using Eq. (1 1), or the estimate using the large po approxima- 
tion [10]. 



y(0) _ r y(0) 


y s,Q 


<ff 


s,2 


vS 


.v.4 


exact (tif = 3) 
Eq.(ll)(n/ = 3) 
large p («/ = 3) 


-1.33333 
-1.23430 
-1.33333 


-1.84512 
-1.95499 
-2.69395 


-7.28304 
-7.53665 
-7.69303 


-37.3395 
-34.0562 


-236.882 


exact (rif = 4) 
Eq.(ll)(n/ = 4) 
large p (»/ = 4) 


-1.33333 
-1.29036 
-1.33333 


-1.64557 
-1.69672 
-2.49440 


-5.94978 
-6.07826 
-6.59553 


-27.6301 

-27.0349 


-161.155 


exact (rif = 5) 
Eq.(ll)( W/ = 5) 
large p («/ = 5) 


-1.33333 
-1.41383 
-1.33333 


-1.44602 
-1.42799 

-2.29485 


-4.70095 
-4.72881 
-5.58246 


-19.4623 

-21.0518 


-103.190 



the detectors correspond to the fields in the Lagrangian (as in QED). Unfortunately, this 
is not the case in QCD and one reflection of this fact is that the pole mass suffers from 
renormalon singularities. Moreover, these renormalon singularities lie close together 
to the origin and perturbative calculations have gone very far for systems with heavy 
quarks. At the practical level, this has reflected in the worsening of the perturbative 
expansion in processes where the pole mass was used as an expansion parameter. It is 
then natural to try to define a new expansion parameter replacing the pole mass but 
still being an adequate definition for threshold problems. This idea is not new and has 
already been pursued in the literature, where several definitions have arisen [11]. We can 
not resist the tentation of trying our own definition. We believe that, having a different 
systematics than the other definitions, it could further help to estimate the errors in the 
more recent determinations of the MS quark mass. Our definition, as the definitions 
above, try to cancel the bad perturbative behavior associated to the renormalon. On 
the other hand, we would like to understand this problem within an effective field 
theory perspective. From this point of view what one is seeing is that the coefficients 
multiplying the (small) expansion parameters in the effective theory calculation are 
not of natural size (of O(l)). The natural answer to this problem is that we are not 
properly separating scales in our effective theory and some effects from small scales are 
incorporated in the matching coefficients. These small scales are dynamically generated 
in n-loop calculations (n being large) and are of 0{me~ n ) (we are having in mind a 
large |3o evaluation) producing the bad (renormalon associated) perturbative behavior. 
In order to overcome this problem, we may think of doing the Borel transform. In that 
case, the renormalon singularities correspond to the non-analytic terms in 1 — 2u. These 
terms also exist in the effective theory. Therefore, our procedure will be to subtract the 
pure renormalon contribution in the new mass definition, which we will call renormalon 
subtracted (RS) mass, mR$. We define the Borel transform of mRs as follows 

B[m RS }=B[mos}-N m v f — -— TTE (l+c i (l-2u)+c 2 (l-2u) 2 + ---), (15) 
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where Vf could be understood as a factorization scale between QCD and NRQCD (or 
HQET) and, at this stage, should be smaller than m. The expression for mRs reads 

r \ ^ a. (fioY n +\t \ x^ T(n+l+b-k) 

^ = mos " 5> v ' (s J a ' (V/) g c ' r(i +t -t) • (16) 

where cq = 1. We expect that with this renormalon free definition the coefficients 
multiplying the expansion parameters in the effective theory calculation will have a 
natural size and also the coefficients multiplying the powers of a s in the perturbative 
expansion relating mRs with m^- . Therefore, we do not loose accuracy if we first obtain 
mRs and later on we use the perturbative relation between m^s and my^ in order to 
obtain the latter. Nevertheless, since we will work order by order in a s in the relation 
between m^s and niy^, it is important to expand everything in terms of a s , in particular 
OL s (Vf), in order to achieve the renormalon cancellation order by order in oc v . Then, the 
perturbative expansion in terms of the MS mass reads 

CO 

«Rs(V/)=m M s+£rf<+ 1 , (17) 

n=0 

where rjp = r^ s (m^g,v, V/). These rjp are the ones expected to be of natural size (or 
at least not to be artificially enlarged by the first IR renormalon). 

In Ref. [1], we have applied this scheme to potential NRQCD and HQET. For the 
former, by using the T(15) mass, we have obtained a determination of the MS bottom 
quark mass. For the latter, we have obtained a value of the charm mass by using 
the difference between the D and B meson mass. In both cases the convergence is 
significantly improved if compared with the analogous OS evaluations. 
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